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We report detailed magnetic, transport and thermodynamic properties of metastable Rri3Bii4 
single crystals in superconducting and normal state. We show that Rli3Bii4 is nearly isotropic, 
weak to intermediately coupled BCS superconductor, whereas the electronic resistivity above super- 
conducting T c — 2.94 K is dominated by the phonon scattering in the large unit cell with pores filled 
by Bi atoms. Superconductivity is strongly influenced by the nature of atoms that fill the voids in 
the crystal structure. 

PACS numbers: 74.25.-q, 74.70.Ad, 74.62. Bf 



I. INTRODUCTION 

Complex crystal structures with atomic subunits of- 
ten show significant electronic and magnetic tunabil- 
ity through incorporation of various elements in the 
(sub)structure. Examples include covalently bonded IV 
group thermoelectric clathrates that are under influence 
of anharmonic lattice vibrations with electron - phonon 
coupling strength comparable to MgP>2, alkali doped 
fullercncs and magnetic molecules*^— Cluster units may 
substantially define physical properties whereas their size 
may considerably deviate from the length scales defined 
by the lattice parameters, even in intermetallic systems 
where delocalized metallic bonds are present £ 

Binary Rh-Bi intermetallic phase diagram has been 
studied for more than eight decades. Yet, structurally 
complex metastable Rh3Bii4 with 136 atoms per unit 
cell has been discovered only in the recent years.— ~— 
Bismuth-rich side of this phase diagram includes also a- 
RhBi 2 , /3-RhBi 2 , RhBi 3 and RhBi 4 phasesJ^ Accord- 
ing to the Rh-Bi phase diagram, Rh.3Bii4 is a metastable 
phase hidden in the RhBi 4 -Bi solid phase region. Since 
the kinetic of peritectic reaction forming RhBi4 is ex- 
tremely slow, it can be easily suppressed by undercool- 
ing due to nucleation problemsAi Hence, Rli3Bii4 single 
crystals can be decanted from the Bi liquid whereas the 
RhBi4 single crystals arc difficult to obtain. 

Rli3Bii4 has orthorhombic symmetry (Fddd space 
group) with a = 0.69041(4) [0.68959(15)] nm, b = 
1.73816(9) [1.7379(3)] nm, and c = 3.1752(2) [3.1758(6)] 
m m 11 ' 12 In a Rli3Bii4 unit cell, there are 24 Rh and 
112 Bi atoms Each Rh atom is surrounded by eight 
Bi atoms, however there are two inequivalent Rh po- 
sitions in the unit cell. Rhl is cubically coordinated, 
whereas Rh2 exhibits a square antiprismatic coordination 
as shown in Fig. 1(a). Each [RhBi 8 / 2 ] polyhedra shares 
common edge with four other polyhedra, and thus create 
a porous three-dimensional (3D) framework with chan- 
nels along the [100] direction. In these channels, there 
are chiral inner surfaces and filled with bismuth ions (Fig. 
1(b)). Partial oxidation of RhBi4 using bromine results 
in Rh3Bii 2 Br2, isostructural with Rli3Bii4 J£ The differ- 
ence between RhjBi]^ and Rli3Bii2Br2 is in the atoms 



in the channels. In Rh3Bii4, Bi atoms fill the channels, 
whereas they are filled by Br atoms in Rli3Bii2Br2. The 
theoretical calculation indicates that the 3D polyhedral 
frameworks in Rli3Bii4 and RhsBi^B^ arc formed by 
covalently bonded rhodium and bismuth atoms, consis- 
tent with the results of high pressure studies^ On the 
other hand, there are no direct bonds between this 3D 
polycation framework and anions (Bi -1 or Br" 1 in the 
channels), which leads to the higher compressibility along 
a axis when compared to those along b and c axes mainly 
controlled by the covalently bonded polycationJ^ 

Magnetization measurement indicates that Rli3Bii4 
is a superconductor with T c = 2.82(5) K^i Although 
the crystal structure and chemical bonding have been 
studied ) 11 ' 12 ' 15 the studies on physical properties, es- 
pecially in connection with superconductivity are still 
lackingAi In this work, we report detailed analysis of 
Rli3Bii4 physical properties in the superconducting and 
normal states. We show that electronic scattering in the 
normal state is dominated by phonons, whereas in its su- 
perconducting state Rli3Bii4 is an isotropic, weakly to 
intermediately coupled BCS superconductor. 



II. EXPERIMENTAL 

Single crystals of Rli3Bii4 were grown by the flux- 
growth method with Rh:Bi = 5:95 molar ratio. Rh pieces 
(99.9 %) and Bi shot (99.9 %) were weighed, combined 
into alumina crucible, covered with quartz wool and 
sealed into the quartz tube with partial pressure of argon. 
The quartz tube was heated to 1000 °C held constant for 
2 h and then cooled at a rate of -3.5 °C/h to 300 °C where 
crystals were decanted. Single crystals with typical size 
0.5x0.2x0.2 mm 3 were obtained. Crystal structure and 
phase purity were examined by powder and single crystal 
X-ray diffraction pattern (XRD) with Cu K Q radiation 
(A = 0.15418 nm) using a Rigaku Miniflex X-ray ma- 
chine. The structure parameters are extracted by fitting 
the XRD spectra using the Rietica software^ Crystal 
was oriented using Bruker SMART APEX II single crys- 
tal X-ray diffractomcter. The composition of Rli3Bii4 
single crystal was determined by examination of multi- 
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FIG. 1. (a) Crystal structure of Rri3Bii4. Three-dimensional 
framework is built of [RhBi 8 / 2 ] cubes and square antiprisms 
which share common edges. Bi4 atoms fill the channels of 
the framework, (b) Edge-sharing [RhBi 8/ / 2 ] cubes and square 
antiprisms enclosing bismuth ions form the zigzag channel in 
Rh3Bii4. (c) Powder XRD pattern of Rli3Bii4 single crystal, 
(d) The EDX spectrum of a single crystal. 



pie points on the crystals using energy dispersive X-ray 
spectroscopy (EDX) in an JEOL JSM-6500 scanning elec- 
tron microscope. Electrical transport, heat capacity, and 
magnetization measurements were performed in a Quan- 
tum Design PPMS-9 and MPMS XL 5. 
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FIG. 2. (a) Temperature dependence of the resistivity p a (T) 
of Rri3Bii4. Inset: resistivity near T c . (b) Temperature de- 
pendence of dc magnetic susceptibility of RhsBii4 with ZFC 
and FC at 10 Oe for H ||a, H ||b and H ||c.The supercon- 
ducting transition temperature T Cl<mse t are marked by arrow. 



III. RESULTS AND DISCUSSION 



According to the phcnomenological model 
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All powder X-ray reflections (Fig. 1(c)) can be indexed 
in the Fddd space group. The refined lattice parameters 
are a = 0.6891(1) nm, b = 1.7388(1) nm, c = 3.1718(2) 
nm with R p = 5.220, R^ = 6.720 and x 2 = 0.919, con- 
sistent with the values reported in literature ! 1 ^ 12 On the 
other hand, we also refined the lattice parameters us- 
ing single crystal XRD. The lattice parameters are a = 
0.707(3) nm, b = 1.755(5) nm and c = 3.154(8) nm, which 
are also close to the result obtained from powder XRD 
and previous results. EDX spectrum of a single crystal 
(Fig. 1(d)) confirms the presence of only Rh and Bi with 
the average atomic ratios Rh:Bi = 3.0(4):13.8(4). 

The temperature dependent electrical resistivity p(T) 
of Rh3Bii4 is shown in Fig. 2(a). The sharp resistivity 
drop with T Ct0nset — 3.04 K is caused by superconduct- 
ing transition (inset of Fig. 2(a)). The residual resistiv- 
ity ratio(RRR), defined as p(295K)/ p(3MK), is about 
9.3. The curve p(T) is convex above 50 K, with a ten- 
dency to saturate at high temperature. The saturation of 
p(T) could be related to the Ioffe-Regel limit ji£ when the 
charge carrier mean free path is comparable to the inter- 
atomic spacing and/or to the two-band conductivity^ 
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where Vf is the Fermi velocity, VL p is the plasma fre- 
quency, a is interatomic spacing and p ideal = Pr + Pi(T). 
The Pi{T) is the inelastic resistivity due to the phonon 
contribution whereas p sa t is the saturation resistivity. 
Phonon contribution can be explained by the Bloch- 
Griinciscn formula: 



ji rQn/T 



(e x - 1)(1 - e~ x ) 



^dx (2) 



where 0jj is the Debye temperature, c is a constant 
which depends on the electronic structure of the metal 
through the Fermi velocity vp and the density of states 
at the Fermi energy.— The fitting result over the full 
temperature range of pr is shown in Fig. 2(a) as a solid 
line, revealing an excellent agreement between the model 
and data for p r = 60.0(4) pfl-cm, p sat = 640(2) pU-cm, 
B D = 95.0(8) K, c— 27.7(6)xl0 4 pQ-cm, 

To confirm the presence of bulk superconductivity in 
Rli3Bii4 single crystals, the magnetization is measured 
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FIG. 3. (a) Magnetization hysteresis loops of RhsBii4 at T 
= 1.8 K for H ||a, H ||b and H ||c. (b ), (c) and (d) Low 
field parts of M(H) at various temperature for H ||a, H\\b and 
H\\c with demagnetization correct, respectively. The solid 
blue lines are the "Meissner line" as discussed in the text, 
(e) Temperature dependence of H c \ for H\\&, H\\h and //|| c. 
The dashed lines are the fitted lines using H c \ = _ff c i(0)(l- 
(T/Tc) 2 ).(f) the temperature dependence of anisotropy of 
ffci, Th 1 = H cl , a (T)/H cl , b (T), 7h , = H cl , c (T)/H cl , b (T), 
7 Hcl = H cl , c (T)/H cl , a (T). 



using dc susceptibility method. Fig. 2(b) shows the tem- 
perature dependence of the dc susceptibility for a crystal 
with dimension of 0.2x0.4x0.2 mm 3 for H \\ a, H || b and 
H || c with zero- field cooling (ZFC) and field cooling (FC) 
under an applied magnetic field of 10 Oe. The steep tran- 
sition show that the Rli3Bii4 single crystal has a super- 
conducting transition temperature for all field directions 
of 2.91 K(T C}0nset ) with a transition width AT C = 0.30 
K, consistent with the result of resistivity measurement. 
At 2 K, the ZFC dc susceptibility approaches a value of 
-1 for all directions after demagnetization correction, in- 
dicating bulk superconductivity. The FC signal shows a 
flux exclusion of 61 %, 70 % and 81 % for H || a, H \\ b 
and H || c respectively The big superconducting volume 
fraction for FC indicates the rather weak flux pinning 
effects. 

The shape of M(H) curve at 1.8 K (Fig. 3(a)) confirms 
weak pinning, in agreement with Fig. 2(b). In order to 
obtain the lower critical fields of Rti3f3ii4, we measured 
the initial M(H) curves at various temperatures with 
the field directions along the a-, b- and c-axis. For each 



field direction, all curves clearly fall on the same line and 
deviate from linearity for different temperature. Linear 
fits for the initial parts of magnetization curves describe 
the Meissner shielding effects ("Meissner line"). 

The value of at which the field starts to pene- 
trate into the sample can be determined by examining 
the point of deviation from the Meissner line on the ini- 
tial slope of the magnetization curve. The first pene- 
tration field H* x is not the same as the real lower criti- 
cal field H c x, due to the geometric effect. The H c i can 
be deduced from the H* ly assuming that the magneti- 
zation M = —H c i when the first vortex enters into the 
sample. Thus H has been rescaled to H = H a — NM 
and H c i = H* x /{1 — N) where TV is the demagnetiza- 
tion factor and H a is the external fields We estimate 
demagnetization factors 0.47, 0.23 and 0.59 for H || a, 
H || b and H \\ c by using H cl = H* cl /tanh(- v /0.366/a), 
where a and b are width and thickness of a platc-likc 
superconductor^ The corrected data are plotted in Fig. 
3(b), (c) and (d). Considering the demagnetization fac- 
tors, the obtained slopes of the linear fitting at the 
lowest temperature of our measurements T = 1.8 K 
are -0.972(2), -0.988(2) and -0.987(2), very close to - 
1 {A-kM = -H) for all field directions. Thus the full 
Meissner shielding effect in our measurement provides a 
reliable way to determine the value of H c \ . 

H c i is determined as the point deviating from linearity 
based on the criterion AM = (M m - M th ) = 2 x 10 -6 
emu, M m is the measured moment value and Mth is fitted 
moment value at the same external field. From Fig. 3(b), 
(c) and (d), H£ lt H h cl and values at the different 
temperatures can be obtained as shown in Fig. 3(e). The 
H cl (T) can be fitted according to H cl (T) = H cl (0)[l - 
(T/T c ) 2 }. Fitting lines are shown in Fig. 3(e). The 
obtained #£(0), # c \(0) and ff^(0) are 90(1), 63(1) and 
93(1) Oe respectively. The anisotropy of H c \ along each 
axis can be determined as shown in Fig. 3(f). It can 
be seen that the anisotropy of H C \{T) is very close to 1, 
especially for H \\ c and H || a, indicating that Rh 3 Bii 4 
is a nearly isotropic superconductor. 

With increasing magnetic fields (Fig. 4(a) and (b)) 
the superconducting transition width broadens and the 
onsets of transition shift to lower temperatures gradually 
for both H \\b and H \\ c. When H = 6 kOe, the super- 
conducting transition cannot be observed above 1.9 K for 
both field directions. It is interesting to observe kinks in 
the R(T) curves below T c under fields for both field di- 
rections. The kinks are becoming more pronounced with 
the increasing field. The origin of these kinks implies 
rich vortex physics, similar to NbSe2.— The upper crit- 
ical field H C 2 is determined by the criterion of 90 % of 
the normal state resistivity at various fields for both field 
directions (Fig. 4(c)). The slopes dH c2 /dT for H || b and 
H || c are equal to -6.47(5) kOe/K and -6.01(11) kOe/K, 
respectively. For dirty limit superconductors, H C 2(0) can 
be obtained from the Wcrthammcr-Hclfand-Hohcnbcrg 
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FIG. 4. (a) Temperature dependence of the resistivity in a 
set of magnetic fields from to 6 kOe for (a) H|| b and (b) 
H 1 1 c. (c) Temperature dependence of the upper critical field 
H c2 for H|| b and (b) H|| c. 



formula^ 



H c2 (0) 



-0.693- 



dH r 



dT 



(3) 



yielding H b c2 {Q) = 13.4(1) kOe and H c c2 {Q) = 12.4(2) kOe. 
Since the Pauli limiting field H p (0) = 18.4T C - 52 kOe^ 
the orbital effect should be the dominant pair-breaking 
mechanism. The small difference in the H c2 values for 
both field directions indicates that Rh3Bii4 shows almost 
isotropic H C 2(T) as seen in Fig. 4(c). 

Since H%y « H^, according to the anisotropic 
Ginzburg-Landau (GL) theory, we assume £ a ~ £ c 
and coherence length £(0) can be estimated from the 
H c2 (0) with: H b c2 (0) = $ /[2^a(0)? c (0)] and H c c2 (0) = 
$ /[27rCa (0)6(0)], where $ = 2.07 x 10~ 15 Wb. Based 
on the values of H c i(0) and H c2 (0), GL parameters 
Ki(0) is obtained from the equation H l c2 (0)/H l cl {0) = 
2K 2 (0)/lnKj(0), where i denotes the field applied along 
i direction. The thermodynamic critical field H c (0) 
can be obtained from H c (0) = F^(0)/[\/2«i(0)].The 
GL penetration lengths are evaluated by the equations 
k 6 (0) = A a (0)/&(0) and n a (0) « k c {0) = A o (0)/&(0) = 
[A a (0)A b (0)/^(0)6(0)] 1/2 . The anisotropy is lams = 
H b 2 {0)/H^ 2 (0) = 6(0)/^(0)^ All of the obtained pa- 
rameters arc listed in Table 1. 




•H 

X 




■ 
# 
■ 
■ 




■ 
■ 
■ 

- ■ 


2.89K 


■ 


. 2.94K 




2.0 2.5 3.0 3.5 4.0 4.5 




T(K) 



10 15 20 25 30 

T 2 (K 2 ) 



FIG. 5. (a) Low temperature specific heat C p /T vs. T 2 from 
1.95 K to 300 K at H = kOe. The solid line shows the 
fit (see text)(b) Temperature dependence of the electronic 
specific heat of RhsBii4 plotted as C es /T vs T at H = kOe. 



Specific heat divided by temperature (Fig. 5(a)) shows 
a jump at 2.94 K, indicating the bulk superconducting 
transition. Specific heat includes both the electron and 
lattice parts C{T) = C e (T) + C ph (T). In the normal 
state, the specific heat of the lattice part is expressed 
by the /3T 3 term at temperatures far below the Dcbyc 
temperature Qjj and electronic specific heat is assumed 
to be - jT. Using C p /T = 7 + /3T 2 we obtain (Fig. 
5(a)) 7 = 8(1) mJ moV 1 K~ 2 and <d D - 113.9(7) K 
using <d D = (12 7 r 4 iVi?/5^) 1 / 3 , where N = 17 is the 
number of atoms per formula unit and R is the gas 
constant. Qd is similar to the value obtained from 
p{T). The electron-phonon coupling constant \ e - p h is 
obtained from the McMillan equation: 



* e— ph 



fi*ln(lA5T c /Q D ) - 1.04 
1.04 + Zn(1.45r c /e D )(l - 0.62 M * 



(4) 



and assuming the common value for the Coulomb pseu- 
dopotcntial fx* w 0.13. The value of X e - P h is determined 
to be 0.74(1) by using T c = 2.94 K and 6^ = 113.9(7) 
K. The value of X e - P h implies intermediately or weakly 
coupled BCS superconductivity. The electronic specific 
heat C es in the superconducting state (Fig. 5(b)) is ob- 
tained by subtracting the lattice contribution estimated 
from the total specific heat. The extracted specific heat 
jump at T c (AC/jT c = 1.8(2)) is somewhat larger than 
the weak coupling value 1.43, and also points to interme- 
diate coupling strength^ 

Rri3Bii2Br2 has similar structure to Rh3Bii4, but 
the two Bi4 atoms are replaced by the Br atoms. 
Rfi3Bii2Br2 is a metal above 2 K without superconduct- 
ing transition^ From theoretical calculation^ the den- 
sity of state (DOS) at Fermi energy is similar for both 
compounds, however the contributions of atoms (s and 
p states of Bi) at Bi4 position to DOS are nearly re- 
moved when voids in the crystal structure arc filled with 
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TABLE I. Superconducting parameters of RhgBii4. 



Rh 3 Bi 14 T c H' cl 


Hl 2 {0) Ht(0) Ki(0) fc(0) 


Ai(0) 7anis 


In 


P 




(K) (Oe) 


(kOe) (kOe) (nm) 


(nm) 


(mJ/mol K 2 ) 


(mJ/mol K 4 ) 




i = a 2.94 90(1) 




273(4) 1.08(3) 


8(1) 


22.39(4) 


0.74(1) 113.9(7) 1.8(2) 


i = b 63(1) 


13.4(1) 5.4(1) 17.4(2) 16.9(2) 


166(8) 








i = c 93(1) 


12.4(2) 6.7(1) 13.1(2) 15.71(6) 











Br atoms. Therefore, the superconductivity in RhsBiu 
might be related to the s and p states of Bi atom in the 
interstitial Bi4 position of Rh3Bii4. 

IV. CONCLUSION 

In summary, we present a comprehensive study of 
the normal and superconducting state properties of 
metastable intermetallic superconductor, RhjBi^. The 
temperature dependence of resistivity p(T) above the su- 
perconducting transition can be explained by the phonon 
scattering. Rli3Bii4 is an isotropic, weakly to intermedi- 



ately coupled BCS superconductor. Presence (absence) 
of superconductivity in Rli3Bii4 (Rh3Bii2Bi"2) suggests 
that the nature of atoms that fill the pores in the struc- 
ture is rather important. 
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